Profile scaling in decay of nanostructures 
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The flattening of a crystal cone below its roughening transition is studied by means of a step flow 
model. Numerical and analytical analyses show that the height profile, h(r,t), obeys the scaling 
scenario dh/dr = T{rt~ 1/4 -). The scaling function is flat at radii r < R(t) ~ t 1/A . We find a 
one parameter family of solutions for the scaling function, and propose a selection criterion for the 
unique solution the system reaches. 
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In recent years it has become technologically possi- 
ble to design and manufacture crystalline nanostructures, 
which are of tremendous importance for the fabrication of 
electronic devices. In many cases, these nanostructures 
are thermodynamically unstable, and tend to decay with 
time. This phenomenon has triggered experimental and 
theoretical efforts to try and understand the decay pro- 
cess ft|-|l3|l. Under fairly robust conditions, the decay of 
a nanostructure at low temperatures (below the rough- 
ening temperature, Tr) is dominated by the motion of 
atomic steps on the surface. Hence, attempts have been 
made to understand and predict the relaxation dynamics 
of simple step configurations. 

In this work we analyze, numerically as well as analyt- 
ically, the time evolution of a crystalline cone formed out 
of circular concentric steps. The decay of other types of 
nanostructures such as bi-periodic surface modulations 
has been studied experimentally on Si(OOl) by Tanaka et 
al. . Rettori and Villain M studied this problem theo- 
retically in the case of small amplitude modulations. Our 
study, on the other hand, is relevant to large amplitude 
modulations and in this sense is complimentary to their 
work. We find that the height of the cone decays with 
time as h(0) — h(t) ~ t 1 / 4 and the radius of the plateau 
at the top of the cone grows with time as R(t) ~ i 1 / 4 . 

Consider the surface of an infinite crystalline cone, 
made out of circular concentric steps of radii ri(t), sep- 
arated by flat terraces. The step index i grows in the 
direction away from the center of the cone. We assume 
no deposition of any new material, no evaporation and 
no transport of atoms through the bulk. To calculate the 
time dependence of the radii, we have to solve the diffu- 
sion equation for adatoms on the terraces with boundary 
conditions at the step edges, taking into account the re- 
pulsive interactions (of the form G(ri+\ — ?"i)~ 2 ) between 
steps. Using a standard approach to do this BR, we 
arrived at a set of equations of motion for the step radii. 
It is convenient to present these equations in terms of 
dimensionless radii, p.;, and dimensionless time r: 
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C® q is the equilibrium concentration of a straight isolated 
step, T is the temperature, T is the step line tension, fl is 
the atomic area of the solid and D s is the diffusion con- 
stant of adatoms on the terraces, k is a kinetic coefficient 
associated with attachment and detachment of adatoms 
to and from steps. 

The equations of motion in terms of these variables 
take the form 
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where the velocities of the first and second steps are mod- 
ified to include interactions only with existing steps. Eq. 
(|l|) depends on two parameters: g and q. g = mea- 
sures the strength of step-step interactions G relative to 
the line tension T, while the parameter q = (1 + jp 1 ^)^ 1 
specifies the rate limiting process in the system. When 
D s 3> k (or q — ► 1), diffusion across terraces is fast and 
the rate limiting process is attachment and detachment 
of adatoms to and from steps. When D s <C k (or q — ► 0), 
the steps act as perfect sinks and the rate limiting process 
is diffusion across terraces. 

We integrated Eqs. (|l|) numerically both in the dif- 
fusion limited (DL) and in the attachment/detachment 
limited (ADL) cases. When the repulsive interactions 
between steps are weak (i.e. g is small), there is a strik- 
ing difference between the dynamics in the two limits. 
In the ADL case the system becomes unstable towards 
step bunching, whereas in the DL case there is no such 
instability. However, when g is large enough the insta- 
bility disappears even in the ADL case. We limit our 
discussion to situations where the step bunching insta- 
bility does not occur. Fig. [l] shows the time evolution of 
initially uniformly spaced steps with unit step separation 
in the ADL case. 
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FIG. 1. Time evolution of the step radii in the ADL case 
with g = 0.01. The front can be fitted by a r 1//4 power law 
(dashed line). 

Each line in the figure describes the radius of one step 
as a function of time. We note that the innermost step 
shrinks while the other steps expand by absorbing the 
atoms emitted by the first step. When the innermost 
step disappears, the next step starts shrinking and so 
on. The disappearance time of the n th step, r n , grows 
with n as r„ ~ n 4 . This process results in a propagat- 
ing front which leaves a growing plateau behind it. At 
large times, the (dimensionless) position of this front be- 
haves as pfront(T~) ~ r 1 / 4 . A similar picture arises in 
the DL case with differences in the details of the indi- 
vidual step trajectories. This power law is an indication 
of a much more general and interesting phenomenon. It 
turns out that for large times, not only the front posi- 
tion but also the positions of minimal and maximal step 
densities scale as t 1 / 4 . In fact, the step density D(p,r) 
obeys the following scaling scenario: There exist scaling 
exponents a, (3 and 7 which define the scaled variables 
x = pT~P/' 1 and 8 = r 1 / 7 . In terms of these variables 
D(p,r) — 8 a F(x,8), where the scaling function F is a 
periodic function of 8 with some period 8q. Our ansatz 
is somewhat weaker than standard scaling hypotheses, 
which would assume F is independent of 8. The necessity 
to introduce a periodic dependence of F on 8 is a man- 
ifestation of the discrete nature of the steps. Thus the 
disappearance time of step n is r„ = (ji8q) 7 . An immedi- 
ate consequence of the scaling ansatz is that if we define 
8 = 8 + n8 with < 8 < 9 , and plot &- a D (xd 13 ^^) 
against x, all the data with different values of n and the 
same value of 8 collapse onto a single curve. 

To verify that our system obeys this scaling ansatz, we 
define the step density at a discrete set of points in the 
middle of the terraces: 
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Pi+i(r) - pi(r) 



(2) 



Fig. |^ is a plot of D{p, t) as a function of x = pr^ 1 ^ for 
a fixed value of 8 and 12 different values of n in the ADL 
case. The excellent data collapse shows that our scaling 
ansatz indeed holds with a = 0, (3 = 1 and 7 = 4. Data 
collapse of similar quality is achieved in the DL case with 
the same values of the exponents. The dependence of the 
scaling function on 8 is very weak in the DL case, and is 
more pronounced in the ADL case. 
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FIG. 2. Data collapse of the density function in the ADL 
case with g = 0.01. The values of the scaling exponents used 
here are a = 0, /3 = 1 and 7 = 4. This figure shows density 
functions with 12 different values of n and the same value of 
8, as a function of x = pT -1//4 . Some of the unsealed data are 
shown in the inset. 

The above results suggest that the time evolution of 
the system can be described by a density, which is a 
continuous function of both position and time. In the 
remainder of this paper we will derive such a continuum 
model, carry out a scaling analysis to evaluate the scal- 
ing exponents analytically, and calculate the scaling func- 
tion. 

Motivated by the simulation results we assume that 
the scaling ansatz holds. This is already sufficient to 
calculate the values of a and (3. First, we derive a rela- 
tion between these scaling exponents by considering the 
height profile h(p,r). Assuming steps of unit height, the 
profile is related to the step density by 



h(p.T) = ho(r) 



D(p',T)dp' , 



(3) 



where ho(r) is the height at the origin. Far enough 
(when p — > 00), h(p,r) does not change with time, i.e. 
lim /3 _K X3 (/i(p, t) — h(fi, 0)) = 0. Combining this with Eq. 
(0), we arrive at the expression 
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h o (0)-h o ( T ) = 



>+r> 



(F(oo,0) - F(x,9))dx , (4) 



where we have used the definition of the function F. On 
the other hand, ho(0) — h (T n ) = n because r n is the dis- 
appearance time of the n th step. 7 satisfies the relation 
T n ~ n 1 , and therefore we have ho(0) — h (T n ) ~ t^ 7 . 
This and the t^ 13 ^ 1 dependence in Eq. (||) lead to the 
relation a + /3 = 1 . 

In addition, conservation of the total volume of the 
system implies that J ph(p, r)dp is independent of r. 
Integration by parts of the derivative of this integral with 
respect to r yields the following equation: 
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Evaluation of this integral in terms of the function F and 
the scaled variables x and 9 shows the integral diverges 
unless a = |pl| . Combining this with a + f3 = 1, we 
conclude that (3 = 1. 

To evaluate the scaling exponent 7 and the scaling 
function F, we continue with the equation for the full 
time derivative of the step density D: 

dD _ dD dD dp 
dr dr dp dr 



(6) 



Eq. (g) can be evaluated in the middle of the terrace 
between two steps (i.e. at p = (pi(r) + p i+ i(r))/2). The 
l.h.s. of @ is calculated by taking the time derivative of 
Eq. (§): dD/dr = -D 2 {p l+1 - Pi ). This together with 
the fact that dp/dr = ((p\{t) + /j i+ i(r))/2 leads to the 
relation 



dD p i+1 + pi dD 



dp 



— +D 2 (p i+1 - Pi ) = 



(7) 



where the step velocities pi can be expressed in terms of 
the pi's using Eq. (0). 

Now we change variables to 9 and Xi = p%9 1 , and 
transform Eq. (|7[) into an equation for the scaling func- 
tion F. In terms of these variables, Eq. @) takes the 
form 



X%-\-\ Xi — 



F((x i+ i + Xi)/2,9) 



(8) 



According to this, the difference between successive Xi's 
is of order 9~ x wherever F does not vanish. In the large 
9 (long time) limit these differences become vanishingly 
small. This allows us to go to a continuum limit in the 
variable x, by expanding all the terms in the equation 
for the scaling function F in the small parameter 
The final result of these manipulations is the following 
differential equation for F: 



~F'-^+ 9^ (f 1 ■ - H 
°-4 + O(*-, = 0. 



F 2 B 



(9) 



A and B are known expressions involving F,F' ,F" ,F"' , 
F"" , where the primes denote partial derivatives with 
respect to x. The existence of derivatives up to fourth 
order in this equation is a consequence of the fact that 
each step "interacts" with four other steps (two on each 
side) through the equations of motion (|l]). A detailed 
derivation of the continuum model and the exact expres- 
sions for A and B will be given elsewhere jl5| . 

Consider Eq. (^) at large 9. Our expansion in the small 
parameter O^ 1 is valid only at values of x where F does 
not diverge or vanish (see above). Therefore, the first 
term in Eq. (^|) is O(l). This term has to be canceled 
by the second term if we require F to satisfy a single 
differential equation. Hence, we must have 



7 = 4 



(10) 



The fourth term vanishes as 9 — > 00, since 7 — 5 < and 
the third term must vanish as well. Therefore, in the 
large 9 limit, F is only a function of x, and we are left 
with an ordinary differential equation for F: 
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F B = 



(11) 



Let us emphasize that our continuum model is valid for 
arbitrarily large surface curvature and slope (unlike other 
treatments [§,^|). Moreover, since our model is an ex- 
pansion in the truly small parameter 9~ x (see Eq. (]§)) 
it becomes exact in the large 9 (long time) limit. Note 
that in going to the continuum limit we lost the periodic 
dependence of F on 8, which is a manifestation of the 
discrete nature of the steps. 

We now turn to study the solutions of Eq. ([n]) . We will 
consider only the DL case, but an equivalent treatment 
can be applied to more general situations Jll| . In the DL 
case Eq. (ll|) becomes [@ 
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g I YIF'F'" + 3FF"" + 9F 
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x~ 2 + ~~)~~~x^~ ( ] 

Without loss of generality, we choose the boundary con- 
dition at infinity to be F(oo) = 1. Any other choice is 
equivalent to our choice with a different value of the in- 
teraction parameter g, since Eq. ( |l2| ) is invariant under 
the transformation F — > aF, x — > a^l^x and g — > aT 2 g. 
Numerical solutions of Eq. ( p"2| ) which satisfy F(oo) = 1 
indicate that there exists a point x — xq near which F 
behaves as F(x) ~ *Jx — xq. An analytical expansion of 
Eq. ([[2]) for small F also leads to the same conclusion. 
Thus, our model naturally predicts a singular point in 
the profile at which F — > 0. We can prove jl5| that xq is 
the scaled position of the boundary of the plateau at the 
top of the hill. 
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In the derivation of Eq. ( pL 2[ ) we assumed F(x) ^= 0. 
This is clearly violated on the plateau. Eq. ( |l2| ) is there- 
fore valid only at x > xq, while for x < xo, the solution 
for F is simply F — 0. Now we have to solve Eq. ( O ) 
for x > xq with some boundary conditions at x = xq 
together with the condition F(oo) = 1. In addition, we 
have to make sure that all the atoms expelled by the 
growing plateau at x < xq are absorbed by the steps at 
x > xq. This is taken care of by the conservation law (||). 
We rewrite (||) in terms of scaled variables and get the 
following equation: 







x z (F{x) -l)dx= x z (F(x) -l)dx- 



J x 

To carry out the last integral we integrate Eq. (jl^) mul- 
tiplied by x 2 from xq to oo, and obtain the boundary 
condition 



lim <^ g (6FF' + xFF" + xF r 
-3x 2 FF"' - 9x 2 F'F") - 



= 



(13) 
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FIG. 3. Data collapse of the density function in the diffu- 
sion limited case with g = 0.01. Simulation data from several 
collapse periods are shown as dots, while the dashed line is 
the xo ~ Xq solution of Eq. (|l2|). 

So far, we have not determined the value of xq. In 
fact, we solved Eq. ( |l2| ) numerically and found a fam- 
ily of solutions satisfying the boundary conditions, which 
differ in the value of xq. For xq < Xq « 1.44 the equa- 
tion does not have a solution. However, for any value of 
> Xq there is a single solution that satisfies the bound- 
ary conditions. Despite the existence of many solutions, 
our simulations indicate that the system reaches a unique 
scaling solution independent of initial conditions. Fig. ^ 
shows an impressive agreement between the xq = Xq so- 
lution and the data collapse of density functions taken 
from the simulations in the DL case. Thus, the system 
dynamically selects the scaling state with the minimal 
value of xq. The precise nature of this selection mech- 
anism is not yet understood and will be investigated in 
the future. 



In summary, we have presented a complete description 
of the relaxation process of an infinite crystalline cone 
below it's roughening transition. The hypothesis that in 
the long time limit the step density exhibits scaling leads 
to an accurate continuous model for the morphological 
evolution of the crystal. Using the model, we were able 
to derive the exact scaling exponents and the differential 
equation that describes the scaling function. This equa- 
tion admits many solutions, and the system dynamically 
selects the one with the smallest plateau. We hope this 
work will motivate new experiments in which our predic- 
tions will be tested. 
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